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The  problem  of  noise  radiation  from  helicopter  rotors  has 
gained  prominence  due  to  its  annoyance  to  the  public  and  detect- 
ability. Alth'^jgh  the  rotor  is  one  of  the  several  noise  generating 
sources  of  helicopters,  it  is  the  most  important  in  the  external 
regions  of  the  present  machines.  Clearly,  the  reliable  prediction  of 
this  noise  in  the  design  stage  of  the  rotor  is  an  important  step  in 
controlling  the  level  of  the  noise  intensity.  There  has  been  a 
steady  advance  in  the  last  decade  in  the  prediction  of  rotor  noise 
(ref.  1).  There  are  still  disagreements  between  the  theoretical  and 
experimental  results  of  rotor  acoustics.  In  addition  to  this  short- 
coming, the  available  theories  suffer  from  a combination  of  the 
following  restrictions: 


Compactness  of  the  acoustic  sources 
Hovering  helicopter 
Observer  in  the  far  field 
Limited  airfoil  shapes 

Limited  surface  pressure  distribution  models 
Singularities  in  the  solution  for  high  rotor  tip  speeds 
Neglect  of  the  thickness  noise 


a. 

b. 

c. 

d. 

e. 

f. 

g- 

It  is  believed  that  the  removal  of  these  restrictions  and  the 
inclusion  of  the  nonlinear  propagation  effects  should  result  in 
reliable  prediction  of  the  rotor  noise. 

Traditionally,  rotor  noise  has  been  divided  into  several 
categories  such  as  rotational,  vortex  and  thickness  noise.  These 
can  be  grouped  into  two  broad  classes  — those  depending  on  the  local 
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pressure  and  viscous  stress  distribution  on  the  rotor  blades  and  those 
due  to  t^ie  nomicl  velocity  distribution  on  the  blades.  For  example, 
rotational  noise  belongs  to  the  first  class  and  thickness  noise  to 
the  second..  A theory  which  incorporates  the  effects  of  surface 
pressure  and  normal  velocity  distribution  on  a moving  body  is 
developed  in  reference  2.  The  formulation  is  then  specialized  for 
propellers  and  helicopter  rotors.  In  this  work  a study  of  compactness 
assumption  of  sources  on  moving  bodies  has  revvealed  that  in  the  case 
of  helicopter  rotors  and  propellers,  the  sources  on  the  blades  cann  t 
be  considered  compact  for  the  observer  position  in  a large  region  oi 
space  around  the  rotor.  If  the  compactness  restriction  is  removed, 
then  one  would  like  to  remove  the  restrictions  of  limited  airfoil 
shapes  and  surface  pressure  distribution  models  to  improve  the 
prediction  technique. 

The  present  paper  discusses  a new  computer  program  developed 
by  the  authors  at  NASA  Langley  Research  Center  based  on  the  results 
of  reference  2.  The  purpose  of  developing  this  program  has  been  to 
remove  the  restrictions  of  the  already  existing  theories  and  thus 
achieve  a new  capability  in  the  prediction  of  the  rotor  and  propeller 
noise.  The  acoustic  computation  is  performed  in  the  time  domain  and 
the  resulting  pressure  signature  is  then  Pouri^r  analyzed  to  get  the 
acoustic  pressure  spectrum. 

Examples  are  presented  in  this  paper  to  demonstrate  the 
capabilities  of  this  new  program.  These  examples  are  selected  mainly 
with  regard  to  the  restrictions  discussed  earlier  which  are  removed 
by  the  new  formulation. 


THE  ACOUSTIC  F0R»’'JLATI0N 

The  formulation  derived  in  re.'crence  2 is  briefly  discussed 
here.  Consider  a moving  body  whose  surface  is  described  by 
■f(y.  t)  = 0 where  t is  the  source  time.  Let  Vp  be  thg  local 
normal  velocity  of  the  surface,  the  acoustic  pressure  p (x,  t)  is 
given  by 
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Where 

X , t: 

c: 


r: 


6: 

P: 


Ti.  Ta*. 


observer  position  and  time 
speed  of  sound  In  undisturbed  medium 
density  of  the  undisturbed  medium 
|x  - y| , y source  location  on  the  body 
the  angle  between  radiation  direction 
r * X - y and  the  outward  normal  to  the  body 
(under  the  integral)  the  surface  pressure  on  the  body 
the  curve  of  the  intersection  of  the  collapsing  sphere 
g = T - t + r/c  * 0 and  the  body  f (y,  t)  * 0 
the  times  when  the  sphere  g « 0 enters  and  leaves 
the  body,  respectively 


For  application  to  rotors  and  propellers,  the  above  equation  will  be 
rewritten  in  the  form  given  beir>w.  Let  a new  frame  n'  be  fixed  to 
each  blade  such  that  ni  nr-plane  contains  the  rotor  disk  and 
ri2-axis  is  along  the  span  of  the  blade.  Let  n?  * T (nt,  nH  and 
11?  = h (ni»  h;)  be  the  equations  of  the  thickness  distribution  and 
camber  surface,  respectively.  The  components^of  unit  radiation 
vector  (x  - ^)/r  and  the  vehicle  velocity  V in  this  rotating 
frame  will  be  denoted  by  (PI,  PI,  PI)  and  (VI,  VI,  VI), 
respectively.  Equation  (1)  can  be  written  as  follows  (ref.  2): 

p (x,  t)  » ^ [Ii  + I2  + la]  + U + Is  (3) 


The  expressions  for  h to  Is  are 
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(8) 


The  symbols  used  in  the  above  expressions  have  the  following 


meaning: 

Dp: 

disk  plane 

T*.  Tj: 

Vi: 

-vt  + 

V$: 

-Vj  - ntn 

fi: 

D: 

rotor  angular  velocity 

2 2 ^21/2 

[1  - rj  + Tl  (1  - r{  )] 

Ap: 

local  pressure  differental  producing  the  lift 
distribution 

Pt= 

pressure  distribution  on  the  blade  due  to  . 
thickness  distribution  alone 

®h  = 

the  angle  between  the  upward  normal  to  the  camber 
surface  and  the  radiation  direction 

Note  that  in  equations  (4)  to  (8),  the  integrations  are  carried  out 
once  along  the  arc  of  intersection  of  the  collapsing  sphere  g ■ 0 
and  the  projection  of  the  blade  planfonns  in  the  disk  plane. 
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COMPUTATIONAL  METHOD 

* 

Equations  (4)  to  (6)  are  evaluated  on  a computer  using  a 
double  nun»erical  integration  followed  by  numerical  smoothing  and 
differentiation  where  required.  Each  of  the  five  terms  are  integrated 
separately.  The  first  three  are  subsequently  differentiated  and  the 
resulting  five  pressure  contributions  are  added  to  obtain  the  pressure 
signature  and  spectrum. 

At  source  t » a sphere  is  constructed  with  its  center  at 
the  observer  location.  Its  radius  R^  is  selected  such  that  its 
circle  of  intersection,  C'»  in  the  plane  of  the  rotor  is  tangent  to 
the  rotor  disk.  From  this  initial  geometry  the  initial  observer 
time,  ti,  is  calculated  from  t-j  * + Ri/c  where  c is  the  speed 

of  sound  in  the  medium.  The  sphere  is  allowed  to  collapse  by  an 
amount  cAt,  where  t is  the  emission  or  source  time.  During  this 
period,  the  helicopter  rotor  is  allowed  to  translate  and  rotate. 

The  resulting  arc  of  intersection  between  the  rotor  disk  and  the 
new  C'  is  swept  point  by  point  in  a counterclockwise  direction 
until  an  intersection  with  a blade  surface  is  detected  or  until  the 
arc  passes  out  of  the  rotor  disk.  When  a blade  is  encountered,  the 
integrands  of  equations  (4)  to  (6)  are  evaluated  and  subsequently 
the  line  integrals  are  accumulated  point  by  point  using  a trapezoidal 
scheme. 


The  collapsing  process  of  the  sphere  g » 0 is  repeated, 
each  time  yielding  a value  for  the  line  integrals  which  are  accumulat- 
ed for  the  source  time  integration  using  simpson  rule.  This  process 
is  continued  until  it  is  detected  that  the  collapsing  sphere  has 
passed  out  of  the  rotor  disk.  The  integration  is  thus  concluded 
for  the  observer  time  t^  and  the  resulting  integrals  are  saved  for 
further  processing.  Successive  points  are  obtained  in  like  manner. 

To  facilitate  numerical  smoothing  and  differentiation  with 
respect  to  the  observer  time  t,  it  is  required  that  the  t^'s  be 
equally  spaced.  Since  the  relation  between  the  observer  time  t 
and  the  source  time  t is  in  general  nonlinear,  an  iteration  techni- 
que is  used  to  obtain  the  initial  radius  Ri  and  the  corresponding 
source  time  where  the  sphere  g » 0 begins  to  collapse.  The 
smoothing  and  numerical  differentiation  which  is  used  are  presented 
in  reference  3.  It  is  based  on  the  theory  of  finite  Fourier  series 
using  sigma  factors  to  improve  convergence  characteristics  and  to 
reduce  Gibbs  phenomenon.  As  a byproduct  of  this,  the  pressure 
spectrum  of  the  acoustic  signature  is  obtained  quite  easily  using 
intermediate  results  of  the  smoothing  and  differentiation  process. 
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EXAMPLES  DEMONSTRATING  UNIQUE  FEATURES 

« 

The  following  examples  are  selected  with  realistic  data  to 
demonstrate  the  unique  features  of  the  developed  program,  Rectangular 
blade  planform  Is  used  In  all  examples.  This  Is  one  of  the  limitations 
of  the  present  program  which  will  be  removed  In  future. 

In  the  first  two  examples,  the  two-bladed  rotor  system  is  4.58  m 
In  diameter  has  a chord  of  0.356  m.  For  the  first  example,  the 
blade  has  an  NACA  four-digit  airfoil  section  of  12  percent  thickness 
ratio.  The  tip  speed  Is  151.3  m/sec.  The  pressure  distribution  Ap 
corresponding  to  this  tip  speed  was  measured  by  Rabbott  (ref.  4)  for 
various  angles  of  attack.  The  angle  of  attack  here  is  8.5°.  The 
chordwise  pressure  distribution  has  a maximum  at  leading  edge  and  the 
spanwise  loading  has  the  familiar  variation  of  Increasing  towards  tip 
and  reaching  a maximum  at  about  90  percent  of  the  radius.  For  this 
example,  a function  of  two  .arlables  approximating  the  pressure 
aistrlbutlon  in  the  outer  40  percent  of  the  radius  was  first  obtained 
and  was  used  as  an  Input  to  the  program.  The  pressure  Pj  due  to 
the  symmetric  thickness  distribution  was  also  obtained  analytically 
using  the  data  given  In  reference  5 and  corrected  for  compressibility 
effect  by  Prandtl-Glauret  rule.  The  observer  Is  10  m from  the  cente*" 
of  the  rotation  and  45°  above  the  rotor  plane.  The  theoretical  pres- 
sure signature  and  the  pressure  spectrum  are  presented  In  Figure  1. 

The  shape  of  the  pressure  signature  Is  considerably  influenced  by  the 
thickness  noise  even  for  such  a high  observer  elevation.  This  was 
found  to  be  true  for  blades  with  blunt  leading  edge.  In  this  and  all 
the  examples  worked  out  so  far,  the  contribution  of  the  expression 

(see  eq.  (6)),  was  found  to  be  of  the  order  of  10  percent  of 
the  thickness  noise  due  to  aii.  This  Is  expected  on  theoretical 

at 

basis.  The  contributions  of  expressions  U and  I;  are  very  small 
compared  to  the  other  terms  except  very  close  to  the  blades. 

The  second  example  has  rotor  tip  speed  of  259  m/sec  (Tip  Mack 
number  « 0.75).  To  utilize  the  measured  data  of  reference  4,  a 
similarity  rule  Is  applied  to  the  blades  of  the  first  example.  T 
obtain  the  same  pressure  coefficient  cp  as  In  the  above  case,  he 
thickness  ratio  varies  along  span  by  the  following  rule  (ref.  6). 

thickness  ratio  , 0.12 

y/}  - Mj^  \/l  - m/ 

where  M2  * Ozf/c  and  Mi  « Oir/c  where  fli  and  a.  a the 
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angular  velocities  of  the  rotors  of  the  first  and  second  example, 
respectively,  and  f is  the  spanwise  distance  from  the  rotor  center. 
The  angle  of  attack  a in  this  example  also  varies  along  the  span 
as  f ol lows  ‘ 


g ^ 8.5 

\/l  - v'l  - Ml* 

where  a is  in  degrees.  Again  pj  from  reference  5 was  corrected 
for  compressibility  effect.  The  observer  is  10  m from  the  rotot'  center 
and  in  the  rotor  plane.  Figure  2 presents  the  pressure  signature  and 
the  spectrum.  The  signature  is  again  considerably  influenced  by  the 
thickness  noise. 

The  above  two  examples  demonstrate  the  use  of  realistic  pressure 
distributions,  airfoils  with  blunt  leading  edge,  and  blade  twist. 

The  third  and  fourth  examples  demonstrate  that  there  is  no 
limitation  on  tip  Kach  numbers.  In  these  examples  a two-bladed 
rotor  of  10-meter  diameter  and  a chord  of  0.4  m is  used.  The  blade 
length  is  1 m and  a biconvex  wedge  airfoil  section  of  6 percent 
thickness  ratio  is  used.  The  angle  of  attack  Is  2.5®.  The  tip  Mach 
number  is  1.375.  Linearized  two-dimensional  aerodynamic  theory  was 
used  to  calculate  Ap  and  py  which  vary  with  spanwise  location. 
Figure  3 gives  the  pressure  signature  and  spectrum  for  the  observer 
50  m from  rotor  center  and  in  the  rotor  plane.  Figure  4 presents 
pressure  signature  and  spectrum  for  the  observer  50  m from  rotor 
center  but  at  45®  elevation  above  rotor  plane.  The  changes  in  the 
signatures  are  striking  but  expected. 

The  fifth  example  demonstrates  the  forward  flight  capability  of 
the  program.  The  helicopter  speed  is  59.2  m/sec  (115  kts).  The  rotor 
system  is  that  of  HU-IH  which  is  14.64  m in  diameter  and  has  a chord 
of  0.53  m.  The  rotor  rpm  is  324.  The  observer  is  22.9  m from  rotor 
center  and  2®  below  the  rotor  plane.  Due  to  unavailability  of 
reliable  surface  pressure  measurements,  only  the  thickness  noise  is 
presented.  However,  it  was  found  earlier  that  at  high  tip  speeds 
and  in  or  near  the  plane  of  rotation,  thickness  noise  is  dominant 
(ref.  7).  This  conclusion  is  born  out  by  comparing  the  calculated 
pressure  signature,  figure  5,  with  the  measured  signature  in 
refererce  8.  The  peaks  of  the  measured  signature  are  higher  in 
magnituce  but  the  deviation  is  less  than  2 db  which  is  considered 
good  agreement  in  acoustics.  The  exact  effect  of  the  inclusion  of 
the  expression  involving  Ap,  which  is  believed  to  be  important  next 
to  the  thickness  noise,  cannot  be  determined  at  this  stage. 
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/ \ CONCLUSIONS 

* h I i 

' “The  present^aper  discusses  a new  theory  and  a computer  program 
for  realistic  calculation  of  acoustic  pressure  signature  and  spectrum 
of  rotor  and  propeller__npise.^  As  seen  from  the  examples  in  this  paper, 
^jmany  of.  the  common  restrictions  of  already  existing  theories  are 
removedjdsing  the  new  theory  which  is  consistent  with  all  previous 
theories.  Only  deterministic  pressure  fluctuations  may  be  used  in 
the  program  at  this  stage  of  development.  This  will  limit  the 
applicability  of  the  program  to  relatively  high  tip  speeds  where  it 
is  known  that  high  frequency  unsteady  pressure  fluctuations  do  not 
contribute  significantly  to  the  sound  level.  There  are  very  few 
blade  surface  pressure  measurements  and  reliable  acoustic  data  avail- 
able to  test  the  theory  in  full.>  Some  comparison  wUh  experimental 
measurements  has  been  given/1n  reference  7 (using  theoretical 
thickness  noise).  Furtheri^mparison  with  the  measured  acoustic  data 
of  a high-speed  propeller  by  Hubbard  and  Lassiter;  (ref.  Bousing 
limited  aerodynamic  data  in  the  blad.i  tip  region  for  acoustic 
calculations)  has  shown  good  agreement  so  far.  One  important  contri- 
bution of  the  new  theory  is  believed  to  be  the  removal  of  the 
compactness  assumption  which  can  introduce  errors  in  acoustic  computa- 
tions. The  new  capability  will  be  used  to  study  this  effect. — 
Already  it  has  been  found  that  in  most  cases  of  interest  one  only 
needs  to  keep  the  two  expressions  Ii  and  I2,  and  in  some  cases 
one  of  these  two  will  give  a good  estimate  of  the  acoustic  pressure 
of  the  rotor.  More  numerical  examples  and  comparison  with 
experimental  data  are  planned. 
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Figure  1.  Example  1.-  Theoretical  acoustic  pressure  signature  and 
spectrum  of  a hovering  helicopter  rotor  for  an  observer  at  45* 
elevation  above  rotor  plane.  Tip  Mach  number  • 0.44. 
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Figure  3.  Example  3.-  Theoretical  acoustic  pressure  signature  and 
spectrum  of  a hovering  helicopter  rotor  for  an  observer  in  the 
plane  of  rotation.  Tip  Mach  number  = 1.375, 
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Figure  4.  Example  4.-  Theoretical  acoustic  pressure  signature^and 
spectrum  of  a hovering  helicopter  rotor  for  a observer  at  45° 
elevation  above  rotor  plane.  Tip  Mach  number  1.375. 
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Figure  5.  Example  5.-  Theoretical  acoustic  pressure  signature 
(thickness  noise  only)  of  a helicopter  in  forward  flight 
(59.2  m/tec,  115  kts)  for  an  observer  2®  below  the  rotor  plane. 
Advancing  tip  Mach  number  = 0.90. 


